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Abstract. The production rate of negative parity baryons was found to be 
much weaker than that of positive states in RHIC experiments. In the present 
paper we show that this suppression is a simple consequence of the coalescence 
dynamics of hadronization. 
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1. Preludium 

The speculations on the nature of quark gluon plasma (QGP) has a long history. 
It was assumed, that due to the large momenta of quarks and gluons in the plasma 
phase, the interaction between them become negligible small as a consequence of 
the running coupling constant. This idea of non interacting massless quarks and 
gluons in a big bag was used by many authors, e.g. the authors listed in Ref. [^. 

However, it became clear soon, that conditions necessary to create such a 
plasma, as depicted in the cartoon (see Fig. [TJ, cannot be fulfilled in the heavy 
ion collisions. The collision time is too short, the volume is too small, the tem- 
perature is too low to produce this massless QGP. Therefore the investigations 
developed in the direction, that what is the structure of the matter produced in the 
heavy ion reactions. (Unfortunately strongly different structures were also named 
"quark gluon plasma" . Thus it is high time to use different names for the different, 
well defined matter structures.) 

One of the most important qualifiers to characterize the matter is the dominant 
degree of freedom. In the original quark gluon plasma investigations the dominant 
degrees of freedom were the massless quarks and gluons. With the realization, that 
at the hadronization stage these quarks interact strongly, the constituent quarks, 
which have large effective mass, were considered the dominant degrees of freedom. 
This development brought up the idea of coalescence hadronization [ 2 . 
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Background of earliest speculations 



Temp = 800 MeV 

time = 1 hour 
Volume = large 



result: 

Quark Gluon Plazma 

with massless and interactionless 
gluons, quarks and antiquarks 

Fig. 1. Visualization of the cooking of quark-gluon soup. 

Here it is important to emphasize, that two different coalescence hadronization 
models were developed. 

In our model it was assumed, that during hadronization both the quark numbers 
and antiquark numbers are conserved leading to the simple and transparent 
quark counting scenario [El EJ- In the other case one assumes, that new quark 
- antiquark pairs are created during the hadronization, and only the net quark 
numbers are conserved [[3]. 

After these original calculations a large number of new publications, dealing 
with different observation, confirmed the validity of the coalescence model [0 |S1 
El QUI- In the present paper we show that the most recently found suppression of 
negative parity baryons is also the direct consequence of coalescence hadronization. 

2. Introduction 

In the early rehadronization studies the main efforts were concentrated on the pro- 
duction probabilities of the lowest baryon multiplets. The structure of the particles 
belonging to these multiplets were similar: they all belonged to the spherical sym- 
metric / = orbital angular momentum state. Experimentally also these low lying 
angular momentum states were observed. Presently, however, an opposite parity 
state (A(1520)) also has been observed experimentally I12j . 
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Fig. 2. The A(1520)/A ratios measured by NA49 (squares) [El, and STAR (stars) 
[I12|. The horizontal line corresponds to a thermal model [I13j. 

Figure 2 displays a compilation of experimental values of A(1520)/A ratio from 
NA49 [E] and STAR [H2 together with the result of a thermal model from Ref. [ 
I13| . The figure clearly shows the suppression of this ratio in central and semi-central 
AuAu collisions, where quark coalescence is expected. 

In the present paper we repeat our earlier calculations of coalescence of con- 
stituent quarks into baryons but now with the inclusion of opposite parity final 
states. We will demonstrate the effect of the symmetry of the internal wave function 
of the produced hadrons on the transition rates. 

In our model the structure of a single hadronization step is assumed as follows. 
In the initial state we have a diquark (M) interacting with the background quark 
system. Due to this interaction the background quarks form a screening cluster 
(A). An incident strange quark (V) will pass this cluster picking up the diquark, 
forming a new baryon A, which leaves the reaction zone. 

For easier understanding we demonstrate this process with an educational model 
and calculate nuclear cross section for the proton - deuteron pick up reaction: 
p + (A + n) — > d + A [El- Here the "proton" plays the role of the strange quark, 
the "neutron" is the picked up diquark, and the "deuteron" is the final state baryon, 
A. Real deuteron has only s-wave and d-wave wave-function component, the p-wave 
state is missing. Since color forces are much stronger than the realistic nuclear forces 
between real proton and neutron, p-wave baryons exist in the nature. Thus we allow 
the existence of the p-wave deuteron in our educational model. 
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3. Simple quantum mechanical model for coalescence process 

Considering an incident proton with momentum k p in the center of mass of the p 
and (A + n) system, the pick-up cross section can be written as follows [ 14 : 

°p+(;M-n)-x2+^(kp>kd) = — \g p +(A+n)->d+A(kp, kd)| 2 ■ (1) 

Vp 

Here the matrix element of the coalescence reaction is determined by the interaction 
potential V np (r n — r p ) and can be calculated as 

3p+(A + „W + A(kp,k d ) = --^||^(^ d ( rn - rp ). e -^(r„ + r P )/2. 



V np {r n - r p ).* A ,„(£,r n ) • e lk » r <>d£dr n dr p , (2) 



The internal wave function of the produced final particle is noted by <^d,z d , where 
Id is the internal angular momentum of the captured neutron in the ground and 
excited state of "deuteron" . After integration over variable £ one obtains 

5 P+ (A +nW+ A(k p ,k d ) = (rn _ rp ). e -^-(r n +r p )/2 

Zttii J 

■ V np {r n - r p ) • V™(r„) • e ik " r "dr n dr p . (3) 
Here the wave function of the neutron bound to the nucleus A is defined as 

i(r„) - / * *a,„(£, r„)d£ (4) 

Introducing new spatial variables R = (r n +r p )/2, r = (r n — r p ) and the momentum 
difference K = k p — kd, we arrive to the expression: 

5 P +(A+„)^+A(k p ,k d ) = J /^(k p ,r)-G„(K,R + r/2)drrfR .(5) 

In the following we calculate this matrix element in eq.JSJ, where the "deuteron" 
and "neutron" parts are given as 

F d , ld (k p ,r) = ftjr).y„ p (r)e- k ^/ 2 
G„(K,R + r/2) = Vn(R + r/2)e 4K - R (6) 

3.1. Deuteron part 

We shall assume that the interaction potential between the incoming "proton" (di- 
quark) and the picked up "neutron" (quark) has the form: 

v ( \ _ j Vi ns ide for r < a , . 

Vn P [r) -\ V outside forr>a , ( '> 
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together with the assumption of V outS ide — ► oo. The interaction range a is expected 
to be in the order of baryon size. 

The radial wave function of the deuteron will be approximated by the spherical 
Bessel functions as 



Ud,o(r) = n d , ■ j ((r/a) ■ 3.14) 
Ud,i{r) = n d ,i ■ ji((r/a) ■ 4.50) , 

where we used the well known spherical Bessel functions, 

j (z) = sin{z)/z 
ji(z) — sin(z)/ z 2 — cos(z)/z . 



(8) 



(9) 



After the first zero of the Bessel functions, x = 3.14 and x\ = 4.50, we shall 
assume the radial wave function to be identically zero. 
Furthermore, the normalization equations 



f 

Jo 



"d.idMVdr = 1 



(10) 



can be satisfied by introducing nd,o = 0.22525 • a 3 / 2 and rid,i = 0.15327 • a 3 / 2 
normalization factors. 

The complete deuteron wave function can be written as 

4>d,iAr)=Yi d A&,<l>)-Ud,u(r) (11) 
With these notations the " deuteron part" of the matrix element has the form 
F d , ld (k P ,r) = ^(r) • V np (r) ■ e^' 2 (12) 
uO (r) , ul (r) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 
Fig. 3. Radial wave functions of baryon with Id = and l d = 1 states at a = 0.7 fm. 
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3.2. Neutron part 

The neutron wave function, which is assumed to model quark wave function inside 
the deconfined region, will be approximated by a Gaussian: 



V>„(R + r/2) = N -exp 
with normalization factor 

9 

N = 



1 



(R + r/2) 2 



jrV*. fl 3/a ' 

The Taylor expansion of this wave function around r = is written as 



(13) 
(14) 



V>(R + r/2) = V(R)+VV>(R)- (r/2) + A^(R) • (r/2) 2 + ... (15) 
Substituting the Gaussian wave function from eq. (|13f) into eq. I|15l) one obtains 



^(R + r/2) 



R J 



1 + -4 ■ R- (r/2) + - A ■ (2 ■ (x 2 .X 2 + y 2 .Y 2 + z 2 .Z 2 ) - s 2 .r 2 )/A 



.(16) 



Let us insert this expression into eq.JSJ 

fir Md 

Z7m 



i^ d (k p ,r) • N-e 



R^ e 'K R 



1 + -| * R • r + \ * (2 * (x 2 .X 2 + y 2 .r 2 + z 2 .Z 2 ) - s 2 .r 2 ) ) drdR (17) 



3.3. Space integral of the deuteron part 
Let us calculate the following integral: 
M d 



ll (kp) = / F ^ (k *" r)dr 

H d , o (0 r , r ) • u ld (r)V np {r) ■ e~^ r / 2 dQ r r 2 dr (18) 



2tt^ 2 

Inserting the Raighley expansion into ea. l|18(l 

M d 



UK) 



OO I 



Yl d fi(9 r ,<l)r) ■ Ud,i d (r)V np (r) 



J2 zZ i l h{Kr/2)Yi%{0rAr)-Yi >m {6 p ,cl> p )dn r r 2 dr , (19) 



1=0 m=-i 

and using the orthogonality relation 



F;* m (0 r , r ) • Yk,n(O r , 4>r)d£l r = Sl,k-5r> 



(20) 



Quark coalescence into opposite parity 



7 



we arrive to the following expression: 

J Id (k p ) = -»'- • / Ud , d (r)t/„p(r)j id (fc p r/2) • r 2 dr ■ Y ldfi {6 p ,cp p ) . (21) 

Zirn J 

This integral has to be multiplied by the first term of the Taylor expansion of 
the " neutron part" : 

B n (K) = Jn-e-W^ ■Y 0<0 (Sl R )-e iKR dn R R 2 dR 

= N- J e- {R/s)2 j {KR)R 2 dR (22) 

Thus the complete matrix element in first order approximation is given as: 

.9 P +(A + „)^+A(k P ,K) = I ld (k p ) * B„(K) (23) 
Choosing the Z axis in the direction of k p , we have 9 = 0, <j> = Q, and thus 

*iM = -^^^} 1/2 -^J^,iAr)V np m i (k P r/2)-r 2 dr (24) 

Thus the production rate, Ai d depends on the bombarding momentum, k p and 
it is determined as 

A ld (k p )=C-\I ld (k p )\ 2 , (25) 
where C is a constant, independent on l^. 

4. Numerical results 

We calculated the production rates A$ and A\ at interaction ranges a = 0.1, 0.3, 0.9. 
The obtained results are displayed in Figs. 4-9. 

The total transition rate, Ri, is obtained by the integration of Ai d (k p ), weighted 
by the distribution function f(k p ): 

POO 

R ld = / f(k p ) ■ A ld (k p )dk p . (26) 
Jo 

Here f{k p ) is the momentum of the proton relative to the center of mass of the 
neutron and screening cluster system. For the present numerical calculations we 
used a Boltzmann function: 

f(k p ) = Exp[-k p /T] (27) 

(We mention here that for coalescence within parton shower [ I15j a narrower distri- 
bution function should be used.) 
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Fig. 4. Space integrated matrix element 
squared Ao(k) in ea. (|25|l at a = 0.3 fm. 



Fig. 5. Space 
squared A\ (k) 



integrated matrix element 
in ea. (|25|l at a = 0.3 fm. 





Fig. 6. Space integrated matrix element 
squared Ao(k) in eq. at a = 0.3 fm. 



Fig. 7. Space 
squared A\(k) 



integrated matrix element 
in eq.(J2SJ) at a = 0.3 fm. 



A(k) 
0.2 




Fig. 8. Space integrated matrix element 
squared A (k) in eci. (|25|l at a — 0.9 fm. 



Fig. 9. Space 
squared A\ (k) 



integrated matrix element 
in eci. (|25)l at a — 0.9 fm. 



The obtained Rq, R± production rates can be connected to the production 
rates of A(1114) and A(1520), respectively. Assuming a fast hadronization at the 
critical temperature, T c , one can directly compare the obtained Ri/Rq ratios to 
the measured A(1520)/A(1114) ratio. Fig. 10 displays our result for Ri/Rq as 
a function of the interaction range, a (solid line). The shaded area indicates the 
experimental result A(1520)/A(1114) = 0.022 ± 0.01 measured by the STAR at 
RHIC [EI]. The dashed line shows the calculated thermal ratio [EI- From Fig. 10 
one can conclude that the Ri/Rq ratio calculated in our coalescence model using 
the interaction range a = 0.8 — 1.0 fm is consistent with the experimental data. 



5. Summary 

It is an inherent property of the coalescence rehadronization model that the pro- 
duction of the A(1520, J p = | ) baryon is strongly suppressed in comparison to 
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Fig. 10. Suppression factor as a function of interaction range, a. 

the production of A(1114, J p = ± + 

This is due to the fact that in A(1520) the orbital angular momentum of one 
of the constituent quark differs by one unit from that of the corresponding quark 
in A(1114). The strength of the suppression depends on the length of interaction. 

From the above consideration one may conclude, that i) in the STAR Au+Au 
reaction a sort of quark matter was formed, meaning that in this case the dominant 
degrees of freedom are the constituent quarks, and ii) such matter was not formed 
in the p + p reaction or in peripheral Au + Au reactions. 

One has to mention, however, that it is somewhat surprising, that in the SPS 
Pb + Pb reaction at y/s — 17.3 AGeV such a suppression also exists. 

Acknowledgment 

The authors acknowledge useful discussions with T.S. Biro and L.P. Csernai. 
References 

1. J. Rafelski and B. Miiller: Strangeness Production in the Quark- Gluon 
Plasma, Phys. Rev. Lett. 48, 1066 (1982); 




10 



J. Zimanyi and P. Levai 



T.S. Biro and J. Zimanyi: Quarkochemistry in relativistic heavy ion 
collisions, Phys. Lett. B113, 6 (1982); 

B. Mullen The Physics of the Quark-Gluon Plasma. Loot. Notes in Phys., 
225 (1985); 

K. Kajantie and H.I. Miettinen: Temperature measurement of quark-gluon 
plasma formed in high energy nucleus-nucleus collisions, Z. Phys. C9, 341 
(1981); 

M. Gyulassy: Formation of a quark-gluon plasma in nuclear collisions, 
LBL-14512, (1982) 

2. T.S. Biro, P. Levai, and J. Zimanyi, Phys. Lett. B347, 6 (1995); Phys. Rev. 
C59, 1574 (1999). 

3. A. Bialas, Phys. Lett. B442, 449 (1998). 

4. T.S. Biro, T. Csorgo, P. Levai, and J. Zimanyi, Phys. Lett. B472, 243 
(2000). 

5. R.C. Hwa and C.B. Yang, Phys. Rev. C66, 064903 (2002). 

6. A. Bialas, Phys. Lett. B532, 249 (2002). 

7. A. Bialas, Phys. Lett. B579, 31 (2004). 

8. V. Greco, CM. Ko, P. Levai, Phys. Rev. Lett. 90, 202302 (2003); Phys. 
Rev. C 68, 034904 (2003). 

9. R.J. Fries, B. Miiller, C. Nonaka, S.A. Bass, Phys. Rev. Lett. 90, 202303 
(2003); Phys. Rev. C68, 044902 (2003). 

10. D. Molnar, S.A. Voloshin, Phys. Rev. Lett. 91, 092301 (2003); Z.W. Lin, D. 
Molnar, Phys. Rev. C68, 044901 (2003). 

11. V. Fricse for the NA49 Collaboration, Nucl. Phys. A698, 487 (2002). 

12. L. Gaudichet for STAR collaboration, ||nucl-ex/ 0307013. 

13. P. Braun-Munziger, D. Magestro, K. Redlich, J. Stachel, Phys. Lett. B518, 
41 (2001). 

14. Rearrangement collision, Section 34., in L.I. Schiff: Quantum Mechanics, 
Second edition, McGraw Hill, New York, 1955. 

15. R.C. Hwa, C.B. Yang, |nucl-th/0401001| 



